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Abstract 

We solve the difference equation with linear coefficients by the Momentenan- 
satz to obtain explicit formulas for orthogonal polynomials. 
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Introduction 


§1 Since the foundational works of Euler IIE710II and Gauss BGaussll the 
importance of the hypergeometric series 


F{a, b, c,z) 


y {a)n{b)n Z 

ho (^)« ” 


with 


1 







(fl)o = 1 and (fl)„ = fl(fl + 1) • • • (fl + n — 1) 

is well-known, because for example nearly all functions important for physics 
are special cases of it. 

§2 Since Euler it is also known, that the hypergeometric series satisfies the 
following differential equation 

z(l - + [c-{a + b + l)z]^ - flfof - 0 

And since Gauss' work on the same it is familiar, that F can be expressed a 
difference equation of second order with linear coefficients. So in the mentio¬ 
ned paper among many other formulas one finds fhe following: 


0 = {c — 2a—{b — a)z)F{a,b,c,z) + a{l — z)F{a + l,b,c,z) — {c — a)F{a — l,b,c,z) 

which, as Gauss calls it, contiguous relation can be interpreted as difference 
equation of second order in the variable a. If one solves this same equation, 
one would obtain an explicit representation for F. 


§3 And this paper is dedicated to the resolution of similar difference equa¬ 
tions. For this we will use the Momentenansatz, tracing back to Euler (see 
IE123I or 1E594I ), which name comes from the similarity to the moment pro¬ 
blem, solved by Stieltjes. See | Stieltjes| . This method will allow the resoluti¬ 
on of fhe following difference equation 


{ax -h a)f{x + 2) = {^x + b)f{x -h 1) -h {jx + c)f{x) 

which contains the example of the hypergepmetric series as a special case, of 
course. At first we want to explain this method in general and illustrate it 
with well-known examples. 


First section: Explaination of the 
Momentenansatz 

§4 It is propounded to us, to solve the following equation 
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(ax + a) f{x + 2) = {fix + h) f{x + 1) + {jx + c) f{x) 

with arbitrary numbers a, a, fi, b,j,c. Or completely general 

n —1 

{ttnX + a„)f{x + n) = {oLkX + %) f{x + k) 
k=0 

where we want to focus mainly on the case n = 2. The generalisation to 
the case n > 2 is neither difficult nor neccessary for our purpose. So all the 
following will hold only for the case n = 2. 


Remark 


§5 Depending on the nature of the variable x the number of the constants 
and therefore the neccessary initial conditions for an unique solution varies. 
If X G Z, so we need two initial conditions. 

If X G C \ Z, so one needs infintely many conditions. Hence we have for 
holomorphic / 


so especially 


/(x + fl) 


OO 


E 

n=0 


/W(x) 
- - 


/(x + 1) 
/(x + 2) 


~ /W(x) 

h 

n\ 


Therefore we can see 

(ax + a) f{x + 2) = {fix + b) f{x + 1) + {jx + c)/(x) 

as differential equation of infinite order. 

From this we see, that number of conditions is indeed infinite, but countable. 
Furthermore we recongnize, that the equation for x G Z is interpolated this 
way. This will be important later. But for now we start with the definition of 
the Momentenansatz. 
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Definition (Momentenansatz) 

§6 Let a equation of the form from above be propounded, then we call the 
ansatz 


b 

f{x) = I R{ty-^h{t)dt 

a 

a Momentenansatz. Here K(f) and h(t) are two arbitrary functions, which 
will later be determined from the concrete problem. 

Remark 

§7 From the Momentenansatz it follows, that we have to determine the fol¬ 
lowing quantities 


I. The function R{t) 

II. The function //(f) 

III. The boundaries a and h 


For this we introduce the following auxiliary equation: 

Definition (Auxiliary Equation) 

§8 If the equation 

{otx + a) f{x + 2) = {px + b) f{x + 1) + {jx + c)/(x) 

was propunded, we call the equation 


t t t 

(xx + fl) J R^^^{t)h{t)dt = (f^x + b) j R^{t)h{t)dt + {■yx +c) j R^^^(f)//(f)df + R(f)^Q(f) 


the auxiliary equation. 


t 

F{t) = I f{t)dt 
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is the integral function of /. Q(f) is another function to be determined. 

Remark 

§9 Although it might seem absurd, to introduce another function Q(f), we 
will see in one moment, that exactly by this we will be able to solve the 
equation. 

Hence it is possible, to derive a system of differential equations from the 
auxiliary equation, that relates the functions h{t), Q(t) and R(t). Thereafter 
it will also be possible, to determine the boundaries. 

Definition (Resolvent) 

§10 We call conditions following from fhe auxiliary equation the Resolvent 
of fhe Momentenansafz. In fhe nexf paragraphs we want to derive these con¬ 
ditions. 

Derivation of the conditions for Q(f), R(t) and h(t) 

§11 Let us start from our auxiliary equation 

f f f 

[ax + a) J = (f^x + b) j R^{t)h{t)dt + {■yx +c) J R^^^{t)h{t)dt +R{tyQ{t) 

If we differentiate this equation with respect to t and use the Fundamental 
theorem of Calculus, we obtain 


(ax + = {^x + b)Rityh{t) + ( 7 x + c)R^-^ it)h{t) + xR^-\t)Q{t)R'{t) + R^(f)Q'(f) 

And after the division byR(t)^~^ 


{ccx + a)R^{t)h{t) = {fix + b)R{t)h{t) -|- {jx + c)h{t) + xQ{t)R'(t) + R(f)Q'(f) 

If we compare fhe coefficienfs of x on both sides of the equation, we find the 
following sysfem of equations: 
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I. ccR^{t)h{t) = ^Rit)h{t) + ^h{t) + R'{t)Q{t) 

II. aR}{t)h{t) =bR{t)h{t) + ch{t) +Q'(f)K(f) 


which is a system of differential equations for the functions h{t), Q{t) and 

m- 

Corollary 1 

§12 The system of equations is under-determined, because one has only 
two equations for three functions. Hence one take one of the theree functions 
completely ad libitum. One has the three following possibilities: 


I. 

h(t) 

is arbitrary, hence 

m 

and 

Q{t) 

are functions of 

h{t) 

II. 

Q(i) 

is arbitrary, hence 

m 

and 

h{t) 

are functions of 

Q{t) 

III. 

m 

is arbitrary, hence 

h{t) 

and 

Qit) 

are functions of 

m 


Which function is put ad libitum, follows from the particular problem in the 
most cases. 


Corollary 2 

§13 If one chooses h{t) ad libifum, so one has two arbitrary constants, be¬ 
cause R{t) and Q(t) appear in their first derivative. 

If one on the other hand chooses one of the other functions ad libitum, only 
one constant appears. 


Remark 

§14 From the last corollary we see, that for the case x G Z, in which two 
initial conditions have to be given for the uniqueness of the problem, we 
can proceed in two ways. We can either take h(t) ad libitum and solve the 
problem directly. Or we choose R(t) or Q(t) ad libitum, with that restriction, 
that the second condition is met automatically. 
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The case a: G C \ Z is a lot harder and one only obtains a special solution, if 
it is not possible, to satisfy all initial conditions at once. 

Finally we want to explain how we can determine the integral boundaries. 

Definition (Determinant) 

§15 We call the equation 

R"(f)Q(f)=0 

the Determinant for the Foundries of the integral 

t 

Corollary 

§16 So if the Determinant has less than 2 solutions, so obtains no solutions 
with the Momentenansatz. 

If it has exactly 2 solutions, so one finds exactly one solution. 

But if it are n > 2 solutions, so one has ( 2 ) solutions of fhe problem. 

Remark 1 

§17 The definition of fhe Determinant immediately follows from fhe auxili¬ 
ary equation 

f f t 

(ax + fl) J = (f^x + b) j R^{t)h{t)dt + {■yx +c) J R^^^{t)h{t)dt +R{tYQ{t) 

For, if one finds two values a and b of t, so that the Determinant 

R^(f)Q(f)=0 

is satisfied, so one finally finds 

b b b 

(xx + fl) J R^^^{t)h{t)dt = (f^x + b) J R^{t)h{t)dt + {■yx +c) J R^^^{t)h{t)dt + 0 

Cl d Cl 
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and therefore 


b 

f[x) = J R^^^{t)h{t)dt 

a 

as solution of fhe equation 

{ax + a) f{x + 2) = {fix + b) f{x + 1 ) + {jx + c) f{x) 

hence of the propounded problem itself. 

Remark 2 

§18 In certain problems it can be convenient, to put instead of 

t 

f{x) = J R^^^{t)h{t)dt 

more general 

t 

f{x) = J R’^+\t)h{t)dt, 

what becomes clear in the peculiar cases and only shortens the calculations, 
but does not alter the procedure in the solution of the problem. The steps are 
always the following: 

1 . Make the Momentenansatz 

2 . Formulate the auxiliary equation 

3 . Derive the Determinant 

4 . Choose one of fhe 3 functions ad libitum an solve the system of equations 

5 . Solve the Determinant 

6 . Determine the constants from the initial condition 

To 4 . we remark, that the function should be taken in such a way, that can be 
done, so that the Determinant has a solution. 

Now having explained everything, we want to illustrate it with examples. 
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Second Section: Examples 


Example 1: The Gamma Eunction 

§19 Let the following equation be propounded 

r(x + i) = xr(x) with r(i) = 1 

We illustrate the procedure with all steps: 

Step 1 : The Momentenansatz 
We put, that we have 


b 

r(x) = J R^-\t)h{t)dt 

a 

Step 2 : Auxiliary Equation 
We formulate the following equation 

t t 

J R^{t)h{t)dt = X J R^-^{t)h{t)dt + R^{t)Q{t) 

Step 3 : Derivation of fhe Resolvent 

We differentiate the auxiliary equation, divide by and find: 

R{t)h{t) = xh{t)+xR'{t)Q{t)+R{t)Q'{t) 

A comparison of fhe coefficients of the powers of x gives the following sys¬ 
tem: 


I. R{t)hit) =R{t)Q'{t) 

II. 0 = +h{t) +R'{t)Q{t) 

Step 4 : Choosing a function ad libitum and solving the system 
Hence the system is easily solved, we choose 

R(f) = t 

The system then becomes 
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and is solved by 


I. h{t) = Q'(f) 

II. h{t) = -Q{t) 


Q{t) = Ce ^ und hit) — Ce ^ 

Step 5 : Solution of the Determinant 
The Determinant is 


0 = fQ{t) = Cfe-* 
For X > 0 one finds the two solutions 

t = 0 and f = 00 
This leads to the intermediate result 


r(x) = —C J 
0 

Step 6: Determination of the constant of integration 
According to the initial condition it has to be 


Hence 


1 = r(l) = -C J e-^dt = -C • 1 
0 


And our solution then is 


C = -1 


r(x) = j e ^dt 
0 
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Remark 


§20 Hence we just derive the integral representation of r(x), even the con¬ 
dition for the convergence of fhe integrals appeared by itself. 

Further we note, that the choice R{t) = t at the same time gave the most 
simple soultion of the difference equation 

r(x -E 1) = xl{x). 

If one considers R{t) in general, one would find 

00 i 

r(x) = I 
0 

where R{t) would just have to satisfy the following conditions 

R(0) = 0 and R(oo) = oo 

That in some way R{t) = fisa convenient choice, will be seen in the following 
examples. 

It illustrates the always reappearing phenomenon in Physics, that the most 
simple solution for a problem with several solutions is always the physically 
relevant one. Furthermore the Momentenansatz 

b 

F(v) = J h{t)R^^dt 

a 

seems to be a promising ansatz also for other physical problems. 

Example 2: Die Legendre-Polynomials 

§21 The Legendre-Polynomials satisfy the following difference equation 

Pn + i { x ) = (2n -E l ) xPn { x ) - nP„_i(x) 

Note, that the difference equation is one for n, not x. The first two Legendre- 
Polynomials are Pq = 1 arid Pi = x. 

Encouraged by the physical principle, we put: 
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b 

P„{x) = J t"h{t)dt 

a 

so put R{t) = t directly. The auxiliary is 


f f f 

J t"^^h{t)dt = {2n + l)x J Rh{t)dt — n J t’'^^h{t)dt + 


and leads to the Resolvent 


1. h{t){e-2xt + l) = Qit) 
II: h{t){f-xt) =tQ'(f) 

with the solutions 


Q(f) = C\/ fi — 2xt + 1 und 
So the Determinant is 


m = 


c 

-2xt + l 


0 = Q(f)r = -2vt + l 

with the solutions 

t = X + \/ — 1 und t = X — — 

This yields the intermediate solution 


x+V x^ — 1 

Pn(x) = C I 

X-Vx2-1 

and because of Po(x) = 1 we find: 


t" 

Vl - 2xt + f2 


df 


X+Vx2-1 


1 = c 


df 


t-Vx^-l 


Vl - 2vf + f2 


= C 


log(\/f2_2xt + l + f-; 


x+Vx^-l 


-Vx^-l 


= ciog(-i: 


Hence 


= Cin 
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and our general formula. 



Pn { x ) 


x+Vx^—1 



X-Vx^-l 


Vl - 2xt + f2 


df. 


Remark 1 

§22 It was not surprising, that the variable x appeared in the boundaries, 
because the difference equations is for n, so a: is just a constant. It could 
indeed be doubted, if Pi(x) = x, because the recurrence equations for P„{x) 
depends on two initial conditions and we only used Po{x) = 0 explicitly. But 
one easily calculates, that our formula yields the correct value. By our ansatz 
we satisfied one condition automatically. Just note, that we could have taken 
evrey Legendre-Polynomial as the other initial condition and would have 
gotetn the same formual then, which is indeed quite remarkable and was not 
to be exspected. 


Remark 2 


§23 The formula 


Pn { x ) 


x+a/x2-1 



x-Vx^-l 


Vl - 2xt + 


dt 


also directly allows a definition of Pn { x ) for n ^ N; and it would be worth the 
effort to check, if this then agrees with the usual definitions in the literature. 
Note, that all values | with n G Z leads to elliptic integrals. But we do not 
want to consider this in more detail here. 


Example 3: The Hermite-Polynomials 

§24 The Hermite-Polynomials satisfy the following equation 


Hn+i{x) = 2xH„{x) — 2nHn-i{x) with H=(x) = 1, Hi(x) = 2x 
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Ansatz and auxiliary equation are determined in the same way as above and 
lead to the following Resolvent 


with solution 


1 . 

II. 


= \Q{t) 

iQ'(t) 

— 2xt 


Q{t) = und h{t) 

Hence the Determinant is 




0 = 


which is solved by 


t = —ioo und t = ioo 

So we have the following rntermediafe solution 


ioo 

H„{x) = C J and for t = iy 

—ioo 

CO 

— CO 

And because of initial condition it has to hold 


CO 

1 = Ho(x) = iC J e^~'^^dy = 2iC\/ne~^^ 

— CO 

and therefore 


So our formula is 


C = 


2i^/n 
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CO 


H„(x) = 

^ ’ 2-sfni^ 




= ^F(/e-x' 

^in 


where F means the Fourier Transform here. 

In exactly the same way one finds for the Laguerre-Polynomials with the 
difference equation 


(n + l)L„+i(x) = (2n + l-x)L„(x)-nL„+i(x) mit Lo(x) ^ 1, Li{x) = l-x 
the formula 


L„(x) 


1 /■! f"e'n 

Wi^Jo 


with the exponential integral Ei{x). Here again the phenomenon occurs, that 
we only needed one initial condition explicitly, to derive the formula. In ad¬ 
dition this formula does not seem to appear in the literature. It is of great 
interest, how this can be explained mathematically and if this can be genera¬ 
lised to other classes of functions. 


Remark 1 

§25 Now we were able to show the method with some important examples 
and see quickly, that the difference equation for the hypergeometric series 
could be treated in the same way. Without great effort one reaches the well- 
known integral representation 

F{a,b,c,z) = - tz)-^dt 

\ J \ f ^ 

which was proved by Euler. Il£366ll . 

And because the hypergeometric function because of the Sturm-Liouville is 
of greatest importance, these explainations are maybe not without use. 
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§26 In addition is was remarkable, that the ansatz R{t) = t, so the most 
simple possible, always led to the correct result and contained all physically 
relevant properties. What is the interpretation of the general ansatz, needs a 
further investigation. But the foilwing principle is almost more important for 
us. 


Priciple of the simplest solution 

If there are more solutions to a mathematical problem, so the most simple 
non-trivial one is always the one important for physics. 
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